Filling the Polar Data Gap in Sea Ice Concentration Fields Using Partial Differential Equations

The “polar data gap” is a region around the North Pole where satellite orbit inclination and instrument swath for SMMR and SSM/I-SSMIS satellites preclude retrieval of sea ice concentrations. Data providers make the irregularly shaped data gap round by centering a circular “pole hole mask” over the North Pole. The area within the pole hole mask has conventionally been assumed to be ice-covered for the purpose of sea ice extent calculations, but recent conditions around the perimeter of the mask indicate that this assumption may already be invalid. Here we propose an objective, partial differential equation based model for estimating sea ice concentrations within the area of the pole hole mask. In particular, the sea ice concentration field is assumed to satisfy Laplace’s equation with boundary conditions determined by observed sea ice concentrations on the perimeter of the gap region. This type of idealization in the concentration field has already proved to be quite useful in establishing an objective method for measuring the “width” of the marginal ice zone—a highly irregular, annular-shaped region of the ice pack that interacts with the ocean, and typically surrounds the inner core of most densely packed sea ice. Realistic spatial heterogeneity in the idealized concentration field is achieved by adding a spatially autocorrelated stochastic field with temporally varying standard deviation derived from the variability of observations around the mask. To test the model, we examined composite annual cycles of observation-model agreement for three circular regions adjacent to the pole hole mask. The composite annual cycle of observation-model correlation ranged from approximately 0.6 to 0.7, and sea ice concentration mean absolute deviations were of order 10−2 or smaller. The model thus provides a computationally simple approach to solving the increasingly important problem of how to fill the polar data gap. Moreover, this approach based on solving an elliptic partial differential equation with given boundary conditions has sufficient generality to also provide more sophisticated models which could potentially be more accurate than the Laplace equation version. Such generalizations and potential validation opportunities are discussed.